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Abstract 

The entanglement between SU{2) <S> SU(2) internal degrees of freedom of parity and helicity 
for reflected and transmitted waves of Dirac-like particles scattered by a potential step along an 
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of nonrelativistic 2D systems, as for instance those for describing single layer graphene, or even 
single trapped ions with Dirac bi-spinor mathematical structure. 
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I. INTRODUCTION 


The growing interest in mapping controllable physical systems onto the bi-spinor structure 
of the relativistic Dirac equation nm has been recently translated into a useful theoretical 
tool for investigating nonrelativistic two-dimensional (2 D) systems of utmost scientific ap¬ 
peal |5KZ|. For instance, the low-energy excitations of a nonrelativistic electrons in the single 
layer graphene are known for exhibiting a massless Weyl spinor structure often related to 
2D Dirac equation solutions. Likewise, on the frontier of the relativistic quantum mechanics 
connections to the solid state physics, the blackhole-likc properties in Bose-Einstein conden¬ 
sates |8], the simulation of Unruh effect in trapped ions [9], the trembling motion and the 
Klein’s paradox for massive fermions in 2D systems [TO] have all been currently investigated 
in both theoretical and experimental scopes. 

On the other hand, quantum information and special relativity have composed the in¬ 
termediate framework for identifying the entanglement/separability information content of 
Dirac bi-spinors, even if the Lorentz symmetry is assumed HU- A recent assertive result 
[T21 ITT] establishes that the SU(2) 0 SU{2 ) group structure of Dirac bi-spinors are as¬ 
signed to a Dirac Hamiltonian written in terms of the direct product of two-qubit operators, 
He, = aP 0 (p ■ <j( 2 )) + maP 0 Ip , from which the free particle solutions of the Dirac 
equation are given in terms of SU(2) 0 SU{2 ) parity-spin entangled states written as 
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where s = 0 and 1 stand respectively for negative and positive frequencies, and p = pp with 
p as a unitary vector. The spinor described by u(p) in the momentum representation is 
related to the spatial motion of the respective Dirac solution, which is coupled to its spin 
and describes a structureless magnetic dipole moment [[13] . For the first qubit, the kets, 
|+) and |—), are categorized as the mutually orthogonal intrinsic parity eigenstates of the 
spinorial particle, and the corresponding total inner product set as |\l/ s (p, t)| 2 = \u(p)\ 2 is 
supported by a normalization factor given by 
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Such a state vector Dirac solution, T s (p, t) is also used to be written as the sum of two 
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bi-spinor structures in the momentum representation given by 
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for which, a most complete description of its intrinsic quantum correlation properties are 
given in Ref. pL3j. Its corresponding explicit parity-spin entanglement structure can be 
straightforwardly extended to simulate all the above preliminarily mentioned physical sys¬ 
tems since their similar Dirac bi-spinor structures may exhibit some analogous spin-spin 
entanglement, or even additional quantum correlations. 

Quantifying the entanglement becomes more interesting whether one considers that the 
simplest and mostly known scheme for introducing external interaction elements is the step 
potential barrier problem. Even locally, it may affect the dynamical evolution of any Dirac 
bi-spinor structure. As it shall be shown, for SU(2)®SU(2) structures prepared as separable 
helicity eigenstates, the local interaction with a step potential barrier destroys the parity- 
spin separability such that the induced entanglement of scattered states can be exactly 
quantihed. 

Some preliminary issues have already investigated the relativistic linear transmission 
mm and the phenomena of planar diffusion of Dirac particles jT8l [2771 through and above 
potential barriers and steps. In particular, the effects of relative phases between incoming 
and reflected or transmitted amplitude components of diffused waves have been accurately 
quantihed mmm- Once that the Dirac equation solutions can be interpreted as two- 
qubit states - one associated to the intrinsic parity and another to the helicity quantum 
number - the projection of particle’s spin on its momentum direction - the corresponding 
parity-spin entanglement associated to reflected and transmitted waves along the potential 
barrier can be systematically quantihed. 

Our aim is, therefore, to quantify the entanglement between intrinsic parity and helicity 
for rehected and transmitted waves of Dirac particles scattered by a potential step-barrier. 
The impinging particles are assumed to move along a collinear arbitrary direction on the 
x — y plane. Given that incident waves can be prepared as separable quantum superposi¬ 
tions of positive and negative helicity eigenstates, the parity-spin entanglement for rehected 
and transmitted bi-spinors shall be determined through the corresponding entanglement 
quantifier. In case of pure states, it is computed through the von-Neumann (vN) quantum 
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entropy which is driven by the x — y incidence angle, 0, and by the relative phase between 
the incident amplitudes. 

The paper is organized as follows. In Sec. II, the Dirac solutions for the step potential 
barrier problem in diffusion and Klein (paradox) regimes are obtained, with the correspond¬ 
ing amplitudes of probability for reflected and transmitted waves given in terms of incidence 
and transmission angles, 6 and O '. In Sec. Ill, the parity-spin entanglement for reflected 
and transmitted waves described as pure states is computed in terms of the vN quantum 
entropy. The transition between scattering (diffusion and Klein zone |231 ST] ) and barrier 
penetration regimes is identified by the entanglement profile, as well as through the profile 
of the averaged chirality of transmitted and reflected waves. The relativistic limits are ob¬ 
tained, and the possibility of including phase differences between hclicity degrees of freedom 
(DoF) is also investigated. Our concluding remarks are given in Sec. IV, where lessons 
concerning the extension of our results to Dirac-like systems exhibiting spin-spin analogous 
to parity-spin entanglement are drawn. 


II. DIRAC SOLUTIONS FOR THE STEP POTENTIAL BARRIER 

The Dirac Hamiltonian for a free bi-spinor structure physically realized by, for instance, 
an electron, is given by 

H d = — ia ■ V + /3m 2 , (3) 

with the 4x4 matrices a = {a x , a y , a z } and f3 satisfying the algebra relations, 

a k a t + a t a k = 2I5 ki , a/3 + f3a = 0, /3 2 = I. (4) 

In the Pauli-Dirac representation the matrices a and (3 are written in terms of tensor products 
of Pauli matrices acting on two subspaces labeled 1 and 2 as 

a- V = a( 1) ®(V-a (2) ), (3 = a* 1 * <g> jf } , (5) 

such that one could interpret the solutions of Dirac equation as a state of two-qubits carried 
by a massive particle whose dynamical evolution is set up by the Hamiltonian Eq. @ H2I25]. 
Noticing that all the calculations for the step barrier scattering problem (which includes the 
particle’s diffusion and the Klein zone regimes) can be carried out in the momentum space, 
the positive energy plane wave solutions of Dirac equation, | ip±), obtained from 0, can be 
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simply described as eigenstates of the helicity operator h = £ ■ p/p 
±1, given by 


, with eigenvalues 
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exp [-i(Et-p-x)}, (6) 


(7) 


and where the kets |1) and |0) carry the intrinsic parity quantum number, as to have P\p) = 
(— l) p \p), with p — 1 and 0 respectively for odd and even intrinsic parity, and the kets \h±) 
are given in terms of cR ; eigenstates |±) through relation (7). 

Considering that interactions with electromagnetic fields could be included in the Dirac 
equation, the presence of an electrostatic constant potential leads to the non-covariant form 
of the Dirac equation as |23J 

= (Hd + V>, (8) 

from which stationary solutions of = (E-V)'ijj, as well as their corresponding bi-spinor 
structure, could be straightforwardly obtained through the substitutions E —>■ E — V and 
P y/{E — V) 2 — m 2 = q into Eq. (ji|. 

For the potential step, 


V(x,y) = 


0 for x < 0, (region A) 


(9) 


Vo for x > 0, (region B) 

the solution is separated into two spatial regions as depicted in Fig. [lj for x < 0, one has 
the solution composed by incoming and reflected waves, and for x > 0, one has only the 
transmitted wave. 

Relative to the notation, the momenta are set as 


Pi —pcosdx + psmdy, ( 10 ) 

for the incoming wave and, respectively, as 

Pr = — pcosQx + p sin 6y, (11) 

Pt = qcosd'x + qsinQ'y, (12) 

for reflected and transmitted waves, with p 2 = E 2 — m 2 , q 2 = (E — Vo) 2 — m 2 , and 9' 
parameterizing the x — y angle of the transmitted wave. 
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Denoting the amplitudes by 4, for the incoming plane waves with hclicity eigenvalues 
corresponding to h = ±1, and by R± and T±, respectively, for reflected and transmitted 
plane waves also with h — ±1, the bi-spinor structure of the corresponding Dirac plane 
wave solutions can be written in terms of 

I il>i{x,y,t)) = {I + \^ + (E,m,6))+E\^_(E,m,9)))exp[-i(Et-p I -x)\, (13a) 

I ipR(x,y,t)) = (R+\i(> + (E,m,n-0)) + R-\i/>-(E,m,n-6))) 

x exp [—i(Et — pn ■ x)] f (13b) 

I ifa{x,y,t)) = (T + \i/; + (E- V 0 , m,0')) +T_|^-(-E- V o ,m,0'))) 

x exp [~i((E - V 0 )t - p T ■ x)}, (13c) 

with 

|lME,m,fl)>= |l)®(±|+>+e"|->) + \/|^|0)®(l+)±e‘' , |-)) • (14) 

The discontinuity of the potential along the rc-axis implies that q sin O' = p sin 6 such 
that the transmission angle can be implicitly written as the relativistic Snell’s law, 



where p = m/E , v = Vo/E, and 0 < 6 < 7t/ 2. At first glance, one could separate the incident 
energy zones into two pieces, one for v < 1, which corresponds to the particle’s diffusion 
regime, and another for v > 1, which corresponds to the Klein zone regime. However, if one 
considers that the transmitted waves could exhibit oscillatory and evanescent behaviors, it 
would be more convenient to divide the problem into three energy zones depending on the 
values of the parameters p and u. In region B above the barrier, the transmitted bi-spinor 
wave, V’t, is oscillatory if 

1 > v + \J{ 1 — p 2 ) sin 2 0 + p 2 (diffusion zone), 

p < 1 < v — J{ 1 — p 2 ) sin 2 d + p 2 (Klein zone), (16) 

or evanescent if 

|1 — v\ < J {1 — p 2 ) sin 2 # + p 2 (tunclling zone) (17) 

which can also be separated into two pieces, v < 1 and v > 1, as mentioned above. 
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The Klein zone exhibits the Klein paradox [23] in which the reflection probability is larger 
than the incidence probability. Such an excess of particle number is compensated by the 
production of an antiparticle number by the transmitted wave at the step potential PSH23- 
The oscillatory behavior of the transmitted wave for both particle’s diffusion and Klein zone 
regimes occurs for incidence angles such that 

si ,A < (1 ~ b!: = sin% , (18) 

1 — fJL z 


For angles 9 > 9 C , the transmitted wave is suppressed by an evanescent behavior as to have 
9 C analogous to the critical angle of the Snell’s law. 

Once the optical parameters are established, the complete solutions given in terms of 


Dirac bi-spinor components, ipj, i/jr and V’t from Eqs. (13), are resumed by 


*Pa(x, y, t) = y, t ) + 4i r (x, y, t), 


(19) 


for region A, and 




for region B , satisfying the boundary (continuity) condition set as 


Mo,v,t) = Vb(o, y,t). 


( 20 ) 


( 21 ) 


which implies that, 

I + -I_ +R+ -R_ = ^E^ )(T+ -T_), 
ei(«-»')( /+ + / _)_ e -«»+»')( fl+ + fl _ ) = (r + + r_), 

The explicit expressions for R± and T± can be computed through 
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with A written in terms of 


A = 


(/i cos 6* -Hsin#) cos 2 #,- 


(24) 


cos 2 0 c — (1 ± a/ fi 2 cos 2 9 c + sin 2 # c ) (cos 2 # + cos 9 \J sin 2 # c — sin 2 #) 
where the + sign in the denominator corresponds to the results for the diffusion zone, and 


the — sign to the results for to the Klein zone. The results from Eqs. (22a 22d) and (23) 
have been manipulated for (oscillatory and evanescent tunneling) diffusion zone(s). By 
changing 1 — v into v — 1, one straightforwardly obtains the corresponding results for the 
(oscillatory and evanescent) Klein zone(s). For a normalized incoming wave, one notices 
that |/+| 2 + |i_| 2 = 1, and that the probability conservation yields 


|P+| 2 + |P_| 2 + [\T + \ 2 + |T_| 2 ] = 1, 


(25) 


J p,x 


where v q>x and v PjX are the x component of the particle velocities, i. e. 


V a x = Va COS O' = ^ 

q ’ q E - V 


cos O' = cos O' 


'(1 


li- 


(1 - z,); 


Vp. x = Vp COS 


9 = — cos 9 = cos 9 1 — /i 2 

E 


(26) 


For sin 2 d c = 1 there is no potential barrier. In this case, |-R+| 2 = |-R-| 2 = 0 . Likewise, 
when the incidence angle reaches the value of the critical angle, i. e. 9 —y 9 C , the potential 
step region B that drives oscillating waves turns into a barrier penetration energy zone for 
evanescent waves. Fig. [2] depicts the curves for reflection and transmission probabilities, 
respectively, |i?+| 2 + \R-\ 2 and ^ [|2~+| 2 + |T_| 2 ], in both diffusion and Klein zones. The 
results are given in terms of the incidence angle through (sin 9) for three different values of 
sin 2 # c , from which one can notice the transition between diffusing and tunneling behaviors 
for 9 = 9 C . The negative values of [|T + | 2 + |T_| 2 ] in the Klein zone corresponds to 
the production of an antiparticle number in order to compensate the growing behavior of 


|f? + | 2 + \R~\ 2 so as to be consistent with continuity equation from (25). 


III. ENTANGLEMENT BETWEEN HELICITY AND INTRINSIC PARITY 

The Kronecker product structure of solution (J6| allows one to interpret Dirac quantum 
states as describing a composed quantum system having two internal DoF: one associated to 
the intrinsic parity, P, and another associated to the hclicity quantum number, h. The en¬ 
tanglement between these two DoF corresponds to a quantum correlation that measures the 














separability between spin and intrinsic parity subsystems - if such Dirac bi-spinor structure 
is not separable, it is due to the parity-spin (P — h) entanglement [ 28] . 

The density matrix associated to each one of the states a = /, R or T can be written as 


Pc = \i^a){lpc\, 


(27) 


with the bi-spinor structure of /, R and T waves summarized by 


\ip a (x,y,t)) = I Tpa) exp [~i(Et - p a • x)], 


(28) 


such that, 

l^a) = 


\A + K c 


[ |1) (8) (a + \ha !+ ) + a-\h a -)) + fc Q |0) ® (a+\h at+ ) - a_| h a -))], (29) 


where 


and 


1 ±a-p a 
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1 + p 


1 —// — //' 


(30) 


(31) 


with v constrained to values given in terms of sin 0 C and p through Eq. (18). 


As mentioned above, the state p a describes a composed system of two DoF in the com¬ 


posite Hilbert space R = TLp <8) Rh, for which the trace of (27) is computed through 

Tr [p a \ = \a + \ 2 + |a_| 2 , 


(32) 


and from which one notices that Pr and pr are not normalized and satisfy the relation 

Tr\p R ] + ^Tr\p T ] = l. (33) 

Entanglement shall be given in terms of the functionals of the normalized density operator 
as to have 


Pa = 


Pa 


Pa 


Tr [Po 


OL 


-\- \a_ 


(34) 


The probability of measuring odd and even parity projections from the state (34), P al and 
P a> o, are calculated through 


Pci = Trp[|l)(l|p^p] 
Pc, o = Tr P [|0) (0| Pa,p] 


1 + 4 

„2 


K- 


1 + «| : 


(35a) 

(35b) 
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where p^ P = Tift [p^f ] is the density matrix reduced to parity space. Probabilities (35) 


depend only on the parameters p and v as through Eq. (31), and are not affected by the 


incidence angle, 9 , in the barrier diffusion process. By the way, the average parity is 

1 


(P)a = TrpfPp" ] = 






(36) 


with V = Diag{+1,+1, —1, —1} the parity operator. As the previous probabilities, the 
mean parity does not depend on the details of the scattering process, namely on 9 or 9'. 

As Pq is a pure state, it only admits a unique Schmidt decomposition that is given in 
terms of the basis vectors of each Hilbert space ra. Tip, for parity, and ' Hh , for helicity, 
and of the Schmidt coefficients of the pure state decomposition. If it exhibits some level 
of entanglement between P and h DoF, then the reduced density matrix of both of them, 
either p p and ph results into mixed states. One infers that a pure state is entangled if the 
two of its reduced systems are mixed. The entropy of entanglement between two subsystems 
that compose the state p% is therefore defined by S a = S[p a ,p} = S[p a) h\ = — ^AfelogA*, 
[22], where S'fp] is the vN entropy of the corresponding density matrix which is computed 
in terms of the p eigenvalues, A*,. Pure states in the reduced representation set the entropy 
of entanglement equals to zero. It implies that the state p% is separable. Likewise, reduced 
mixed states lead to non-vanishing values to the vN entropy, which implies that p ^ is 
entangled. 

The reduced representation of p a to intrinsic parity subsystem p Qj p = Tr h[_Pa\ has eigen¬ 
values 






KI 2 I 

CK— 
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(Kl 2 + 

| OL— 
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(37) 


The entropy of entanglement is therefore given by 


Sa ^a:,+ lo§2 ^ct,— • (38) 

Fig-1 shows the plots for the vN quantum entropies for reflected and transmitted waves, 
Sp and St, as function of the incidence angle, 9 , in terms of sin 9 , for the same values of 
sin 2 0 c as depicted in Fig. [ 2 } One can notice the discontinuity of the entanglement quantifier 
derivative at sin 9 = sin 9 C , which sets the point for the transition between diffusion and 
barrier penetration zones. It is important to emphasize that, for oscillatory regimes, i.e. 
for 9 < 9 C , the Sp plots return the same profiles, which does not depend on the choice of 
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9 C , and all the corresponding curves are superposed. Once the transition to the evanescent 


occurs, one has 9 > 6 C . By taking a look at the denominator of A from Eq. (24), one 
notices that \/ sin 2 # c — sin 2 0 is then automatically converted into i \J sin 2 0 — sin 2 0 c , i. e. the 
real and imaginary parts of A are changed such in a way that the Sr plots shall return 
different profiles for the barrier penetration zones. Given that such a transformation of 
real(imaginary) into imaginary(real) parts of A is not analytic, it results into a discontinuity 


in the derivatives of the curves obtained from Eqs. (37 38) for the entanglement quantifier. 

Finally, these results can also be extended to the discussion of statistical mixtures as 
Pi = / + | , 0+)('0+| + I- 1'0-)('0-|, for which, the entanglement could be computed either 
in terms of concurrence m or in terms of logarithmic negativity BH. Our preliminary 
calculations have revealed that, for maximally mixed incident states (/+ = /_ = 1/2), the 
entanglement for both incident and reflected states cancels off and vanishes. 


A. Correspondence to the averaged chirality 


For completeness, the chiral profile of the Dirac bi-spinors can be identified through the 
averaged value of the chiral operator, (y 5 ), with y 5 written in the two qubit SU( 2) <g) SU( 2) 
representation as y 5 = cxi 1 ) ® lf\ Explicit calculations yields 

<7 5 }i = (39) 

for the incident wave and 


< 7 % 

(7 5 )t 


±w ,v~ . (in 


(1 — p?) sin 2 6> c + p? 


|r_| 2 ), 


(40) 


for reflected and transmitted waves, with + and — signs in the expression for (y 5 )^ corre¬ 
sponding to diffusion and Klein zones. Fig. [d] shows the behavior of (y 5 )r,t as function of 
sin 9 for the same parameters of Fig. [2] From Fig. [2j one can identify a very nice qualitative 
correspondence between the averaged chirality and the vN quantum entropy from Fig. |3j 
the results for reflected waves in the oscillatory regime do not depend on the choice of the 
critical angle, 9 C , in the sense that the plots are all coincident; the transition to the barrier 
penetration (evanescent) behavior is accurately quantified by 9 = 9 C \ at least qualitatively, 


11 












in both diffusion and Klein zone regimes, the plots indicate the possibility of there being 
some (quantum) correlation between (y 5 )^ and S Rj t- 


B. Extremal points of the entanglement 

In the diffusion and Klein zones, it is possible to analytically compute the extremal points 
that maximize and minimize the entanglement for the reflected states. The extremal points 
satisfy the condition given by 

as R 

= o, 


dd 


e=e 0 


which implies that 
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2 1 

dd 

L(lfl+ 
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= 0. 


(41) 


(42) 


e=e o 


If 1 .1 is real, the above equation has two solutions 


sin^Q 1 * = 0, 
sin^Q 2 ’ = 




y/l + /l 2 


(43) 


The corresponding eigenvalues of p Rj p read 


Ai 11 =\± | Vl - 4(1 - a 2 )|/+ /-I 2 , 
\ (2) _ 1 ± At 

— o > 


(44) 


from which the entropy of entanglement is calculated through (38). In spite of not obtaining 


analytical expressions for the extremal points of the transmitted wave entanglement, one 
can identify the extremal points for the quantum entropy from the plots depicted in Fig. [3j 

C. Relativistic and non-relativistic limits 


The ultra-relativistic limit of the above obtained results are computed through the limit 
of p —>• 0, for which one has 

i sin 6 cos 2 6 r 


A(p —* 0) — Aim — 


cos 2 9 c — (1 ± sin 2 0 c ) ( cos 2 # + cos 6 \/sin 2 6* c — sin 2 0) 
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(46) 
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which indicates that, in the massless limit, the entanglement for reflected and transmitted 
states does not depend anymore on the incidence angle, 9, and equals to the entanglement 
for the incident state. 

Otherwise, the non-relativistic limit given by fi —> 1 leads to 

A(p —> 1) = Amr = e ld . (47) 

as well as 

K R = K T n 

(i + 4) 2 (i + 4) 2 ’ 1 J 

A +t R = 1, \-j{ = 0 and the parity-spin entanglement is null for all incident, reflected and 
transmitted waves. 

D. Introducing phase differences 

Relative phases can be relevant in a large set of quantum mechanical experiments, for 
instance, in the context of optical interference phenomena, such as in the double slit exper¬ 
iments which include a transparent plate before one of the slits [323- Fc» r the coefficients I± 
described by complex amplitudes, one may write 

I± = e^\I ± l (49) 

such that the relative phase of interest is written as Ac o = oj + — cu_. This phase has no 
effect over |i?+| J + |iiL| 2 or over \T + \ 2 + |TL| 2 > however, it changes the values of |i? + | 2 |i?_| 2 
and |T + | 2 |T_| 2 , as one can notice, for instance, from the explicit expansion of the reflection 
amplitudes, as given by, 

|R + | 2 |R_| 2 = (Im[A]Re[A]) 2 + [(Im[A]) 2 - (Re[A]) 2 ] 2 |/ + /„| 2 

—2(hn[A]Re[A])[(Im[v4]) 2 — (Re[A]) 2 ] sinAcn|/ + /_|(|/ + | 2 — \I~\ 2 ) 

—4 (Im [A] Re [A]) 2 sin 2 Aca|/ + /_| 2 , (50) 

Once that |R + | 2 |i?_| 2 explicitly depends upon the relative phase, it obviously affects the 
results for the entanglement computed through the vN entropy. Results for Sr and St for 
three different values of the relative phase, Acn, are shown in Fig. [5j for the same set of 
parameters used in Figs. [2] and [3j 
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In particular, for a maximal superposition between positive and negative hclicities ob¬ 
tained by setting \I + \ = |/_| = 1/V2 (c. f. Fig.§, the minimal point for the quantum 
entropy in the diffusion zone can be analytically computed for any phase difference, as it 
occurs at 


sin da = 




%/TT 




(51) 


i. e. the same minimal point found in (27). 


Finally, all the above results for the quantum entropies computed from Dirac bi-spinor 
structures related to particles (positive intrinsic parity) can be converted into similar results 
for antiparticles (negative intrinsic parity) : one can perform the substitutions of E —>• —E 
and V —V, and equivalently, p —» — p. One notices that it leads to Re [A] —» — Re [A] into 


Eq. (50) (c. f the third term of Eq. (50)), and the entanglement for the antiparticle state 
is the same as that for the particle state, with opposite phase differences. For Acu = 0 the 
results are equalized. 


IV. CONCLUSIONS 

In this work, the entanglement between helicity and intrinsic parity for 2 D projections of 
3 D plane wave solutions of the Dirac equation with bi-spinor structure has been quantified 
in the context of the step potential barrier problem. Results for reflection and transmission 
coefficients in terms of incoming helicity amplitudes, incidence angles, 6 , and critical angles, 
6 C , for the transition between oscillatory and barrier penetrating evanescent regimes have 
been obtained for impinging particles in diffusion and Klein zone energy regimes. For pure 
state configurations of incident, reflected and transmitted waves, the vN quantum entropy 
was set as the quantum correlation quantifier for the parity-spin entanglement, which was 
computed as a function of 6. Extremal points have been obtained in order to show that a bi¬ 
spinor quantum correlation structure can be manipulated between minimum and maximum 
values through the manipulation of the geometry of an optical-like apparatus that eventually 
controls reflected and transmitted waves. Still in the context of the formal framework, a 
very suitable qualitative correspondence between the averaged chirality, (7 5 )r,t, and the 
vN quantum entropy, Sr,t, have been identified in both diffusion and Klein zone regimes. 
Figuring out the quantitative origin of such an eventual correspondence between chiral 
conversion and entanglement, even in an enlarged context involving Dirac equation solutions, 
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may indicate an interesting routine to the subsequent works. 

Turning our conclusions to a phenomenological view, some recent issues have considered 
the trapped-ion physics [33, CM] as a flexible platform to map several suitable effects in rela¬ 
tivistic Dirac quantum mechanics. Therefore, the planar diffusion and the 2D scattering of 
SU (2) <8) SU (2) bi-spinor structures may also become a useful tool for quantifying elemen¬ 
tary quantum correlations between such relative compounding spinor substructures, or even 
between two arbitrary Dirac particle subsystems. In particular, it has been shown that the 
2D relativistic scattering for ^-dependent potentials may exhibit some kind of entanglement 
between transmitted and reflected wave packets and the transverse momentum [33], which 
could also be quantified through the systematic procedure described in this paper. 

Our final conclusion is that the framework presented here can be straightforwardly manip¬ 
ulated to compute spin-spin entanglement of nonrelativistic 2D systems, in particular, either 
for electron-electron or electron-hole pairs in the single layer graphene or, for instance, for 
single trapped-ions supported by a Dirac bi-spinor dynamics, which indeed deserve a deeper 
investigation. 
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FIG. 1: Schematic picture of the 2D Dirac diffusing process. A Dirac bi-spinor collides with a 
potential step at x = 0 with incoming momentum pj making an angle 9 with the step. Reflected 
and transmitted waves have momentum pr and pr, respectively. 
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FIG. 2: Reflection (first plot) and transmission squared amplitudes (second plot) for diffusion 
(thick lines) and Klein (thin lines) energy zones, namely v < 1 and v > 1, as function of sin (9 for 
fj, = 0.5, 1+ = 1 and sin0 c = 1/2 (continuous line), sin0 c = 1 /\/2 (dashed line) and sin e c = V 3/2 
(dot-dashed line). 6 = 9 C indicates the transition between oscillatory and evanescent behaviors and 
corresponds to a saturation point in the evanescent regime, for which \T\ = 0 and \R\ = 1. 




FIG. 3: Quantum entropies (entanglement), Sr (first plot), and St (second plot), as functions 
of sinR Again, the plots are for /r = 0.5, I + = 1 and sin^c = 1/2 (continuous line), sin@ c = 
l/\/2 (dashed line) and sin0 c = y/3/2 (dot-dashed line) for diffusion (thick lines) and Klein (thin 
lines) energy zones. One notices the discontinuity of the entanglement quantifier derivative at 
sin# = sin0 c indicating the transition between diffusion (c. f. oscillatory behavior) and barrier 
penetration (c. f. evanescent behavior) zones. The oscillatory regime shows coincident results for 
reflected waves in both diffusion and Klein zones. The two horizontal dotted lines correspond to 
the minimum and maximum values of the entropy given by Sr. 
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0.5 



FIG. 4: The averaged value of the chiral operator (y 5 ) as function of sin@, for the reflected (first 
plot) and transmitted (second plot) waves for diffusion (thick lines) and Klein (thin lines) energy 
zones, and with sin# c = 1/2 (continuous line), sin# c = 1 /\/2 (dashed line) and sin Oc = v/3/2 
(dot-dashed line), in correspondence to Fig. [ 2 } Again, 6 = 9 C indicates a (discontinuous derivative) 
transition between oscillatory and evanescent behaviors. 
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0.8 



FIG. 5: The vN entropies, Sr (first column) and St (second column) for \I + \ = |/_| = l/y/2, with 
parameters and plot-style in correspondence with Figs. [2] and [3j The relative phases are given by 
Aw = 7t/4 (first row), Aw = n/3 (second row) and Aw = n/2 (third row). One notices that the 
local minimum for St decreases as the phase difference decreases. For Aw = vr/2, 9 = 9q, with 
sin 6*o = (J,/y/l + fJ. 2 , Sr vanishes in the diffusion zone. 
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